BOUNDARY CONDITIONS FOR SCALED RANDOM MATRIX 
ENSEMBLES IN THE BULK OF THE SPECTRUM 



A.V. KITAEV AND N.S. WITTE 

Abstract. A spectral average which generalises the local spacing distribution 
of the eigenvalues of random N X N hermitian matrices in the bulk of their 
spectrum as N — ► oo is known to be a T-function of the fifth Painleve system. 
This r-function, t(s), has generic parameters and is transcendental but is 
characterised by particular boundary conditions about the singular point s = 0, 
which we determine here. When the average reduces to the local spacing 
distribution we find that r-function is of the separatrix, or partially truncated 
type. 



1. Motivations 

Recent studies [El HU [15] have revealed that a generalised spectral average for 
the universality classes of scaled random hermitian matrix ensembles - the bulk, 
hard edge and soft edge classes for /3 = 2 - are determined by general transcen- 
dental solutions of the Painleve equations Py , Pin and Pn respectively. These 
works unified earlier studies [19] , [23] , [24] , [25] in the sense that this average includes 
local eigenvalue spacing distributions and moments of the characteristic polynomi- 
als as special cases, and extended them in that the parameter sets of the Painleve 
equations were exhausted. As a consequence of this identification the logarithmic 
derivatives of the spectral averages are Jimbo-Miwa-Okamoto u-functions satisfy- 
ing second-order second-degree ordinary differential equations. It was noted in the 
first mentioned works, that while identifying the precise parameters of the Painleve 
transcendents involved, the algebraic approach employed could not uniquely specify 
adequately the boundary conditions that the particular solutions to the differential 
equations should satisfy. It is the purpose of this paper to rectify this for the bulk 
scaling case and completely determine the boundary conditions for the differential 
equations of the cr-form. 

The plan of our work is as follows - in Section 2 we review the isomonodromic 
formulation of Pyi for generic values of the parameters, give a specific parameteri- 
sation of its monodromy matrices and the local expansion of the r-function about 
the fixed singularity t = 0. In Section 3 we give the corresponding expansion for a 
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generalised spectral average of the circular unitary ensemble, the spectrum singu- 
larity ensemble, and identify its monodromy data by comparison with the previous 
section. In Section 4 we review the necessary isomonodromic formulation for Py in 
a parallel manner to the treatment of Pyi in Section 2. We take the bulk scaling 
limit of our spectrum singularity average via a formal argument in Section 5 and 
identify the resulting monodromy data applying in this case. We show that this 
data is consistent with the rigorous theory of the limit transition from Pyi to Py as 
developed in Kitaev [20]. We wish to advise the reader that we re- use the same 
symbols in the context of the monodromy theory of Pyi as for that of Py where 
there is no risk of confusion in order not to burden the notation unnecessarily. 
Where both are discussed together then we make a notational distinction. 



2. ISOMONODROMY DEFORMATION FORMULATION FOR P 



VI 



Following the conventions and notations of [18] , [20] we consider the Lax pair of 
linear 2x2 matrix ODEs for ^(A;t) with four regular singularities in the A-planc, 
denoted by v £ {0, t, 1, oo} 



(2.1) 



dX 



A A x A t 
X A-l X-t 



dt 



X-t 



We adopt the convention that the residue matrices A v satisfy 



(2.2) A +A t + A 1 =-A oo =-^a 3 , a 3 



1 

-1 



6 C\Z, 



and the constraints tr/l„ = 0, dctyl„ 



Ifl2 



2 ,, v £ {0, t, 1}, defining the formal 



exponents of monodromy 0„. The T-junction for Pyi is defined as 



(2.3) 

and the a -function as 



dt 



logr = Tr 



Ao M 
t t-l 



At, 



(2.4) C ( t )=t(t-l)-logT + - 



which satisfies the second-order second degree differential equation 



dt 



(2.5) -C [t(t-l) — C 



dt 2 



dt 



-M-o\) 



- r x + -Aet + e oc y\[-c + -(6t- 



dt 



dt 



5 C + I (f, » 
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Figure 1. Monodromy representation of the fundamental group 
for C\{0,M,oo} 

Furthermore we suppose that the matrices A v ,v G {0,i, 1} are diagonalisablc, i.e. 
that there exist invertible matrices i?„ € SL(2,C) such that R~ l A v R v = hO^a^, 
9 V G C\Z. In the neighbourhood of a regular singularity *(A) can be expanded as 

oo 

(2.6) *(A) = J2 *—( A - "r*" 4 ^, 

m=0 

for v G {0, 4,1} and for A = oo in the form 

(2.7) *(A) = (i + * TOOO A- m J A"^ 3 . 

The monodromy matrices M v (y G {0, t, 1, oo}) arc defined by ^\ V+Se 2*i = ^\ v +s M v , 
and are given in terms of the connection matrices C v by 

(2.8) M v = CZ l e" i8 »° s C v , C x = I. 

They satisfy the cyclic relation which in our conventions is taken to be 

(2.9) M oo M 1 M t M = I, 

and corresponds to the particular basis of loops given in FigfT] 

The monodromy data A^vi : ~ \Qv>C vy M. v \v = 0,t, l,oo} are preserved under 
deformations with respect to t. The invariants of the monodromy data are de- 
fined to be p M = 2cos7r# M := TrM M , /i G {0,t, l,oo} and p^ v = 2cos'K(J t _ LU := 
Tr M^Mu, n,v€{0,t,l}. 

Now we are in a position to present the details of a parameterisation of the mon- 
odromy matrices and the expansions of the T-function about the singular points 
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0, 1, oo given in Jimbo's study [IS] . In this work Jimbo states the following condi- 
tions under which his results apply 

(2.10) 0O,0t>01,0aogZ, 

(2.11) < K(aot) < 1, 

(2.12) ±6t±tTot, e oc ±6 1 ±<j ot ^2Z, 

and we call these generic conditions. In our application the non-resonant condition 
(|2.10p will be adhered to, however (|2.12[) can be relaxed in a meaningful way and 
we conjecture that the results of Jimbo still hold under suitable limiting procedures 
where the left-hand sides — > 2Z. 



When cot 7^ the parameterisation of the monodromy matrices employed by Jimbo is 

I e -iri9^ I ' 



(2.13) M 



, 01 ., 1 / cosTra-e-^cosTT^i -2re- 7r,;e - sin f + (9i + cr) sin § (0^ + X -a) 

(z.l4J iwi 



2 ' 

ismnO^ \2 r - 1 e 7rie °° sinf - 0i + cr) sin f - 0i - cr) -cosher + e 7 ™*^ costt0i 



(2 15) DM D- 1 - 1 I e^cosTT^-cosvr^o -2se™ sin §(0 O + t - cr) sin §(0 O - 9 t + cr) 

i sin ira l v 2s- 1 e- 7r?;<T sin § (0 O + t + a) sin | (0 O - t - a) -e- ma cos 7r0 t + cos 7r0 o 

_j if e" la cos 7T0 O - cos 7T0 t 2s sin § (0 O + 0* - cr) sin §(0 O - 0t + v)^ 



(2.16) DAIqD~ = 

where 



ism-ira \ -2s _1 sin |(0 O - t - cr) sin §(0 O + t + cr) -e^ 7 ™"" cos7T0 o + cos7T0 4 



(2.17) D 



sm-^tfoo -V-L-cr) r sm ^ (^t/oo + + crj \ 
v r- x sin 1(0^ - 0i + cr) sin 1(0^ + 0i - cr) )' 

with the short-hand notation a = a ot , s = s ot . The quantity s ot together with <7ot defines a unique solution to the Pyi cr-form 
(|2.5p . The other parameter r is an arbitrary non-zero complex constant and is a free parameter in the monodromy matrices but 
doesn't appear in the local expansions. 
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A key formula is the following connection relation which relates sot, Pot t° "ti 
and col- 

Lemma 2.1 ([TS],[5]). Under the above generic conditions (2.1012.1 112. 12]) and 

notations 

(2.18) 4s ±1 sin ~(0„ + t T tr M ) sin |(0„ - t ± a ot ) 

7T 7T 

x sin-^oo + 6»i =F o-Qt) sin -(0^ - 0i ± cot) 

_ e ±7Mcr 0t S i n7r( j 0t C0S 7TO~ti — COS 7T0 t COS 7T0QO — COS 7T0 O COS 7T0i) 

± i sin 7rcrot cos 7T(7oi + cos 7r0 t cos ~k9\ + cos 7T0OO cos ttOq. 

A consequence of this relation is a constraint on the monodromy invariants 
which is an algebraic variety defining a sub-manifold, the monodromy manifold, of 



Lemma 2.2 ([18|). The monodromy manifold for Pyi is given by 
(2.19) Wl(p ot ,pti,Poi) =PotPtiPoi +Pm+Pti + Poi 

- (PoPt +PlPoo)P0t - {PtPl +P0Poc)Ptl ~ (POPI +PtPoo)P01 

+ Po+Pt +Pi+pL+ PoPtPlPoo -4 = 0. 

Theorem 2.1 ([18]). Under the conditions \2.10\2.11\2.12\ then we have the as- 
ymptotic expansion of the r-function as t — > in the domain {t G C|0 < \t\ < 
s, |arg(i)| < (/>} for all e > and any <p > 



(2.20) r(t) ~ const • ^^-^-e?)^ 



X \ 1+ S^ 2 1 

l2 (a i „\2\\n1 (a i ~\2l 



_ 1 [g 2 -(g t + q) 2 Pg o -(fli+^) 2 ] ,i- CT nr ,, a(1 _ R(ff)) . 

16<t 2 (1 -a) 2 + Ul J 



where a ^ and s are related to s through 
(2.21) s = s 



r 2 (i - <r)r(i + §(0 O + e t + <r))r(i + ±(-0 o + 0* + *)) 



r 2 (l + <r)r(l + §(0 O +6 t - a))T(l + ±(-0 o + t - a)) 

r(i + i(0oo + 0i + tr))r(i + U- e oo + h + CT )) 



x 



r(i + i(^oo + 0i- o-))T(i + i(-0oo + e 1 - a)) ' 

and we employ the short-hand notation s = Sot,s = §ot and cr = erot. The mon- 
odromy data defining a unique solution to the sixth Painleve system is {<7ot,sot}. 
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3. The Spectrum Singularity Ensemble 

A fundamental ensemble in random matrix theory is the ensemble of finite 
rank (rank = N) random unitary matrices - the Dyson circular unitary ensemble 
( CUE). Consider a member of this ensemble U with eigenvalues z\ = e l6 ' , . . . , Zn = 
gifljv _ Then the eigenvalue probability density function for the ensemble is the Haar 
measure for U(N) 

(3.1) pjv((9li ... j<9Ar ): = 1^ J] l^-^l 2 - 

^ ' ' l<j<k<N 

A generalisation of this ensemble, referred to as the spectrum singularity ensemble 
(SSE) 15], has an eigenvalue probability density function containing additional 
algebraic singularities at z = 0, — 1, — 1/t where t = e l ^,cf> £ [0, 2n). In the log- 
gas picture of the CUE the eigenvalues of a random unitary matrix are mobile 
unit charges subject to a pair-wise mutual logarithmic repulsion and in the SSE 
these charges are subject to additional external fields - the logarithmic electrostatic 
potential of impurity charges located at the ends of the sector, at z = —1,—1/t 
with charges u>i , [i respectively, and an external electric field of strength u>2 ■ The 
generalised spectral average we wish to investigate is the partition function of such 
an electrostatic system, and is given by the Af-dimensional integral 15 



(3.2) AN{t;oj\,W2,^\C) ■= 




N 

x []zp- 2 ii +^1^11+^1^ n i*j-**i a . 

{=1 l<j<k<N 

where (* £ C and the parameters W\, u>2, A 4 £ C, lu = io\ + iu>2, are restricted with 
3?(2a;i), 3?(2/x) > —1, N G Z> . The independent variable t, whilst defined on the 
unit circle \t\ = 1, can be analytically continued into the cut complex t-planc. 

It was shown in p~5] that the average AN(t) is the TV £ Z>oth member of a 
sequence of classical r-functions of the sixth Painleve system. Thus this average is 
characterised by solutions of the nonlinear differential equation (|2.5p subject to the 
boundary conditions given in [17] , which were derived following an idea introduced 
in pi]. 

Theorem 3.1. |17j For generic values of the parameters fjb,u>,uj, subject to 2/i + 
2co>i ^ Z. 5R(2/i + 2uji) > the spectral average An has the following expansion 
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about t = 1 

(3.3) ^w - n r(1 + fc + M+w)r(1 + fc + /i + di) |i 

+ ^W (1 -* )+0((1 -* )) 

+ (-l)^ 1 / e -^(i*-a)) ^ S in7r2^sin7r(/i + cj) 
sin7r(2/i + 2wi) \ 2i sin7r(2^i + 2cji) 

r(i + 2 M )r(i + 2^)r(i + M + o,)r(i + ^ + g) _ , 1+2ai+2wi 

T 2 (2^ + 2wi + 2)r(2^i + 2wi + l)r(7V)r(-7V - 2/z - 2wi) 1 ' 

x (l + 0(l-f)) + 0((l-£) 2+4 ^ +4wi )j- 

The precise relationship between the spectrum singularity average -4/v(i) and 
the isomonodromy theory of the sixth Painleve system was first given in p~7] . How- 
ever because we intend to apply the Kitaev theory [20] , which treats the coalescence 
of the regular singularities 0, t of the Pyi system under the transition limit t — ► 0, 
and the natural bulk scaling of A^(t) is through the limit t — ► 1 the monodromy 
data of our application will have to be recast in a suitable form. 

Theorem 3.2. |17j The spectrum singularity ensemble can be defined by the fol- 
lowing monodromy data. The formal monodromy exponents are 

(3.4) 6 Q = N + 2fi, 9 t = -N-2u 1 , 9 1 = -(x-lo, 9^ = /i + u>, 

with the necessary classical condition 9q — 6 t + 9\ — 9 X G 2Z. The relevant mon- 
odromy invariant is 

(3.5) aot = 2a* + 2wi. 

In this case the associated monodromy coefficient is given by Sot = ens/2 under the 
limiting process e := aot + 9\ — 9^ — > 0, along with 

(3.6) s sm 7r2wi sin 7r(/x + cj) = : h £ 



sin7r(2/i + 2u>i) 2i 
All monodromy matrices are upper triangular 

(3.7) M = 

(3.8) M t = 

(3.9) Mi = 



— 7TZtlo 


mo 







,TTi8 t 


m t 





e —iriO t 


— TTlOi 


mi 
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_ e -™(2 M +2 W1 ) sin7r 2/isin7r(/i + uj) 



+ sin7r(2/i + 2uj\) sin7r2tJi sin7r(/i + uj) 



(3.12) 



mi = — 2«re 



where r is a non-zero arbitrary constant. 

Proof. In |17j it was noted that the monodromy exponents could be given by any 
one of three sets 



modulo permutations of the exponents, an even number of sign reversals and subject 
to some constraints which we elaborate on later. It was found that the first set lead 
to upper triangular monodromy matrices, the second to full monodromy matrices 
except for one which was a multiple of the identity and the third to lower triangular 
matrices. The second set strictly violates the condition of non-resonant monodromy 
exponents and the third set is essentially equivalent to the first and therefore we 
choose an example from the first set. These constraints, when applied to the first 
set, imply that our choices are now 



with an even number of sign reversals. The monodromy invariants were given as 
cr t = N - (i + uj, C7 t i = 2/i + 2wi, (Jqx = N — /x + u). 

Next we employ the linear fractional transformation of our original system t >—> 
1 — t. 8 <-> 0i, cot * "til o~ti l— > °oti an d it is this new system to which we shall 
refer to henceforth. We take one choice from the above (|3. 14|) which is applicable 
to the theory given in and this is (|3.4j) and (|3.5p . Having made this choice 
we observe that in order for the Jimbo parameterisation given in (|2.20|) and (|2.21[) 
to be consistent with ()3.3j) under the map t i— > 1 — t then sot must vanish in the 
manner given in Theorem (|3.2p . The resulting finite coefficient s is then given by 
formula (|3.6|) . It can be verified also that the overall prefactors and the leading order 



(3.13) 




(3.14) 



{6*o, 0t,9i, 9oo} — < 



'±(ji + u), ±{N + 2//), ±(JV + 2wi), ±(fx + Q) 
±(p + Q), ±(N + 2wi), ±{N + 2/i), ±(m + w) 
±(M + w), ±(N + 2/i), ±(N + 2wi), ±(/t + tD) 
±(p + w), ±(N + 2wi), ±(JV + 2/i), ±(/i + w) 
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analytic terms in both expansions agree precisely with the choices we have made. 
The monodromy matrices are then given by p.7H3.12|) using the parameterisation 
in (|2. 1412. 1512. 1612. 17|) . The off diagonal elements satisfy the relation e~ wie °mo + 
e -Ti0t m + e -7ri(e +e t -e oo ) mi = ag requ i rec i b y tnc cyc ij c identity (|2~9|) . □ 

4. ISOMONODROMY DEFORMATION FORMULATION FOR Py 

The work of Jimbo [18j has a slightly different formulation of the isomonodromy 
problem for Py from that of Kitaev and collaborators [20l [H [2] and it is the latter 
form that we adopt. We formulate the Py isomonodromic system as the Lax pair 
of linear 2x2 matrix ODEs for "J (A; t) with two regular singular points at v = 0, 1 
and an irregular one at oo with Poincare rank unity 
(4.1) 

d ft A ^i \ T c? fA 1 



We require that the residue matrices satisfy diag(A + A\) = — %^er 3 , 9^ e C\Z, 
and the constraints trA^ = 0, det^4^ = — i6„, v G {0,1}. Again we assume 
that there exists invertible matrices R v € SX(2, C) such that R~ 1 A„R L , = ^6^03, 
6 V e C\Z, v e {0,1}. 

Within the works of Kitaev and collaborators there are two different conventions 
for the monodromy data employed and we will make the distinction between the two 
by using a caret for those of [20] as opposed to those of p] [2] . In the neighbourhood 
of the irregular singularity the canonical solutions ^(A) can be expanded as 

(4.2) * fe (A) ^ocA- m ) exp (-A - log \)a 3 

m=l ^ 

in the sectors — 37r/2 + irk < arg At < tt/2 + nk, k = 1,2. The Stokes matrices Sk 
are defined as ^ fc+1 (A) = ^(A)^, which have the structure 

(4.3) S 2 

where Sk are the Stokes multipliers. The monodromy matrix is given by 

Mkoo = SkSk+ie 7T ' ie °° a:> , Coo = I, and we set M m = M 0oo - The monodromy 
matrices at the regular singularities arc defined in the same way as in (|2.8| . The 
cyclic relation in this case is 

(4.4) Mo Mi Moo = /, 

reversed in order from the usual convention because of the nature of the limiting 
transition we are to consider later. The monodromy data preserved here are My '■= 
{9 V , C v , M„, si, S2}. If we define the monodromy invariant 2cos7tct = Tr(MoMi), 
then the following constraint, analogous to (|2.19p . applies 

(4.5) sohe-™ 9 ™ =4 sin I (61^+ a) sin I (^-o-). 



* 












: 


! 021+1 — 


K 


1 / 


V s 2 i 
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In contrast the cyclic relation adopted in [H [2] is 

(4.6) MooMiMq = I, 
where 

(4.7) Moo = M 2oo = S 2 e™ e ~ a3 S 1 . 

The inter-relationship between the two sets of monodromy data is given by 

(4.8) S = S U S! = S 2 

(4.9) M = SiMiMoM^S^ \ Mi = SiMiSf 1 , = S , iA'/ 00 S"r 1 . 

The (T-function in the Jimbo formulation of the Py linear system is equivalent 
to that of Andreev and Kitaev [2] and is defined as 

(4.10) CW =*^logr+i(fl + floo)t+ \{(0o + 0co) 2 
which satisfies the second-order second degree differential equation 



dt 



d j. 



dt 



dt^ ~^ 



^dt^ 



-I 2 



1 



dt 



The asymptotic expansion of the T-function is given by the following theorem. 

Theorem 4.1 ([H]). Under the conditions O , 0i ^ Z, < 5R(ct) < 1 and 8 1 ±8 ± 
cr, Oqo ± a ^ 2Z £/ien we Ziaue t/ie asymptotic expansion of the T-junction as t — > 



(4.12) r(i) ~ const •t (<r2 - ( >»)/ 4 i 1 



9 2 + ( t 2 )_ 



4cr 2 



. [floo-a][fl 2 -(fli-a) 2 ] 1+g 
8ct 2 (1 + ct) 2 



-i[^oo+a][0 2 



8ct 2 (1-ct) 2 
where a ^ and s are related to s through 



k 1 - + 0(\t\ 2 ^-^) , 



(4.13) 



r 2 (i - a)r(i + + flp + a))r(i + - O + a)) 
'r 2 (i + <7)r(i + + O - <r))r(i + i(0i - O - a)) 



r(i + i(0oo + g)) 
r(i + i(0oo -*))' 



When a ^ Z we have the parameterisation of the Py monodromy matrices employed in [18j . and in Theorem 6.2 of [2] 
(4.14) 

x _ 1 / e^coSTrflo -CO8 7T01 -2rse- 7r4,T sinf (6>i - 6>o - o-)sinf (6»i +6» 

~ isinTrcr \ v 2(rs)~ 1 e" tT sinf(6' 1 -0o + cr)sinf(e'i+6io + (T) -e^'" 7 cos7r(9 + costt^ 

(4 15) UJlf D- 1 - 1 f e^cosTr^ -COSTCO 2rssin f (0 X - O - <r) sin f (0„ + ft - 

1 isinTrcr ^^(rs)" 1 sin § (ft - ft + ct) sin §(ft + ft + cr) -e^ 7 ™ 7 cosTrft + cosTrft 

where 

/ e -W2 rsin +cr)\ 



»— 1 p ivi<j/2 



(4.16) D = 

The Stokes multipliers are given by the formulae 



r(i-2^)r(z^)' r(i-2^)r(2^«)' 

The parameter r is again an arbitrary non-zero complex constant and doesn't appear in the expansion formulae. 
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5. The Bulk Scaling Ensemble 

Also studied in |15j was the scaling limit of the spectrum singularity ensemble 
in the neighbourhood of the singularity at (j> = (t = 1), to the bulk regime via the 
limit N — > oo. It was shown there that the log-derivative of the average -4jv(i) con- 
verged to a solution of the Jimbo-Miwa-Okamoto a-function for the fifth Painleve 
equation. Thus the most general universality class of hermitian random matrix 
ensembles in the bulk scaling limit was found to be a generic four parameter class 
- three arbitrary parameters, u)i,u>2,/jl, appearing in the differential equation and 
one, £*, in the boundary data. Thus this class is characterised by a transcendental 
solution to a generic fifth Painleve equation. 

The establishment of this result was based on a formal scaling limit of the 
Pvi second-order second-degree ordinary differential equation (|2.5[) to the corre- 
sponding Py ODE (|4.1ip . Defining the scaling variables 

(5.1) t := e~ x/N , u(x;uJi,uj2,fi;C) : = lim log.Ajv(<; Wi, w 2 , £t; £*), 

it was found that u{x) is related to a solution of an alternative Jimbo-Miwa- 
Okamoto cr- form of the fifth Painleve equation h v (s) by 

1 3 

(5.2) u(ix;u>i, w 2 , m) = h v (x; v) + (— —w — -u> + fi)x + y 8 (w — a)) 2 — fx(ui + ui), 
with the Okamoto parameters 

(5.3) (vi,v 2 ,v 3 ,v 4 ) = (wi - i«w 2 , -u\ - ^iuj 2 ,H + ^ iuj 2, + ^«*> 2 ). 

The alternative Jimbo-Miwa-Okamoto a-form of the fifth Painleve equation (|4.11[) 
is 

4 

(5.4) {xKf - [h v - xh' v + 2(h' v ) 2 } 2 + 4 JJ^v + O = 0, 

k=l 

where the constraint v\ + v 2 + v 3 + -u 4 = applies. The scaled spectral average is 
a t- function for this system and related to u(x) by 

(5.5) A(x)=exp — u{y;wi,W2,fi). 

Jo V 

We have the following expansion of A(x) about x = under the above formal 
bulk scaling limit. 

Theorem 5.1. Under the conditions Si(2/z), §R(2o;i), + cj), + id) > — 1, 
< 3?(2^i + 2u>\) < 1, the spectral average -4jv(t; u>i, u>2, {J>; £*) has a bulk scaling 
limit with 1 1— > exp(— x/N) as N — > oo 

ckc^ tt 1 T(l + fc + /i + w)r(l + A; + + oi) 

(5 ' 6) U fc!r(2 M + 2^ + fc + l) * N (t; W ,^ ) nZoo A(x), 
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and the scaled A-function has the following expansion as x — > 

(5-7) A(x) = l + !0^-x + O(x*) 

+ 1 ft* e-^-^ + sin7r2^sin7r(/i + cj) 
7T \ 2i sin 7r(2/i + 2cji) 

v r(i + 2 At )r(i + 2^)r(i + » + ^)r(i + M + g) 1+2m+2lu1 r 

r2(2 M + 2 Wl +2)r(2 / x + 2c 1 + l) 35 u + 

+ 0(.T 2+4 " +4Wl ). 

Proof. The expansion (|3.3[) is valid uniformly in N under the substitution t 
exp(— x/N) and by taking the limit N — > oo and employing the asymptotic formula 
for the ratio of gamma functions we arrive at (|5.7p . □ 



Comparison of this result with the general Theorem 14. II of Jimbo leads to the 
following conclusions. 

Theorem 5.2. The Py monodromy parameters for the bulk scaled spectrum singu- 
larity ensemble are 

(5.8) 9 q = li + u), 6 1 = -/jl-lj, 000 = 2(1 — 2u\, 

with the invariant a = 2\i + 2lo\ . The monodromy coefficient is 

, r „. sin7r2/j siri7r(2/x + 2wi) e -Mn-o) 

(5.9) S = : : : ; r£ — . 

sm7r2o;i sm7r2o;i siri7r(/i + u)) 2i 

Proof. The identification of (|5.4|) with the parameter set (|5.3|) with (|4.11|) leads to 
the following solution sets for the Py formal monodromy exponents 

(5.10) {^o + t^oo: -^o + ^^oo, - 4^00, ^t^ 1 - 4 °°J' 

modulo permutations. One choice, consistent with subsequent calculations, is given 
in (|5.8|) . Using this choice and comparing the two asymptotic expansions (|5.7p and 
(I4.12|) through their relationship 

(5.11) A(x) = (7 e 3( e =c-w+<s)x s i[(eo+eoc.) 2 -e2]-i(8 o +i0 M ) 2 -2 AlWl + i( w -a) 2 r ^ :I ,^ 

we see that a 2 = (2/x + 2uii) 2 from consideration of the algebraic prefactor. The 
coefficients of the analytic terms then agree upon using this value for the exponent. 
Finally the coefficient of the non-analytic x 1 ~ a vanishes because 6q — B\ — a = 
and (|5.9p follows from a matching of the remaining term. □ 
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6. The Limit transition Pvi to Pv 

The limit transition from Pyi to Py is from one monodromy preserving system 
to another and it is known [5] that there exists no continuous transition between 
such systems which is itself monodromy preserving. However there are discrete 
transitions which are and one example is a sequence of even parity Schlesinger 
transformations of the Pyi system, i.e. 9 U i— > 8 V + 2n,n 6 Z + . In a sense the 
Py system arises then as a fixed point of this map although the deformation variable 
also must scale in an appropriate way. This approach was fundamental to the study 
of Kitaev [20] who considered an example which is relevant to the present work. 
By way of comparison with our Theorem 1 5. II let us now consider the formal scaling 
undertaken in limit II of Pyi as given in |20j . In this limit the Pyi formal exponents 
of monodromy are related to the Py exponents by 

(6.1) e tvi = 6 6 -2n=:~, 

(6.2) #ovi = #oov — #6 + 2n = 9oov H — > 

e 

(6.3) 9 lvl = 9 lv , 

(6.4) #oovi = #0v, 

and the limit n — > oo or e — > is taken. The quantity 8$ is an additional constant 
which will fixed in our application. The two regular singularities v = 0, t coalesce 
through t vl = ei v = 0(e) and the transition from the linear Pyi system to the 
Py system takes place according to the scheme displayed in Fig. [5] 

At the level of the Pyi isomonodromic system the following scaling takes place 
as n — ► oo, e — » 

(6.5) A VI = -L, 

(6.6) t V i = ei v , 

(6.7) lim R^ vl (t vl , A VI ) = * v (i v , A v ), 

(6.8) Rtvi2 9 °°vi<r3 R tvi = 4>v + 0(e), 

(6.9) OivAv. = A lv + 0(e), 

(6.10) -R^\J^R tvl = j- QfloovCTa + ^ov + A lv j + 0{e). 

Consequently Theorem 2 of Kitaev [20] allows us to compute the Py monodromy 
data under the limit transition II described above. 

Theorem 6.1 ( ED] ) • Let /i VI £ -M V i(0oov — &6, #ivi> &6, #oovi) an d suppose the 
following conditions hold: 

(1) fPp .ftO )) and 8 uvl e C\Z, 

(2) (M tvi - e-^HMovi - e~^ e ~ v ~ ^) ± 0, 
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Poincare index P V i Pv Pv Poincare index 




oo oo i-> 

Figure 2. Coalescence of the two regular singularities 0,t of the 
Pvi system into the irregular singularity at oo of a Py system 
according to the limit transition II of [20] . 



(3) a, (3,1 G C\Z, where 
(6.11) T = Tr(M tvI - e" e6 )(A/ 0vl - e ~^oov-e e )^ 



(6.12) I = -L log ^coSTrfloov + ± ^(cos^v + iT) 2 -! ] , 



(6.13) 



(4) t/ie inverse monodromy problem for H2.1]) is solvable for all pairs (fi V i-,t VI ) 
such that t VI = et v > 0, 

(5) it is possible to construct the sequence A^^etv), n G Z + . 
Then the data /z v G M v (6ov, 0i v , 0oov) are given by the formulae 
(6.14) 6 0v = 0oo VI , 9i v = 6*i VI , 

(6-15) 5o= L X Sr= I' ™ e , 

\r(i-a)r(i-/3) V \ u 1 / 

(6.16) M 0v = KM ODV1 K-\ M lv =KM lvl K-\ 
where K is the unique (up to a sign ) solution of the system 

(6.17) KM tvl K~ x = S e 7Tl9ea3 1 KM^K' 1 = e^ 6 ^ 3 5 1 e 7riflooVO ' 3 . 

Applying the above Theorem to the spectrum singularity ensemble we have the 
following result. 
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Corollary 6.1. In the limit transition of Theorem W.l\) we find the Painleve V 
parameters of the spectrum singularity ensemble are 

(6.18) 6q v = h + u), 1v = —fi - u, 0QOV = 2/U — 2cji, 0$ = — 2u>i. 

The Stokes multipliers are 

7wi o 7ri „7ri(2 M -2^ 1 ) 



GC 



The monodromy matrices are 
(6.20) 

e ^(M+s>)(i _^*) _ e Ti(M-a) 2 j s i n7r 2 At 



M n 



(6.21) 



(6.22) 



r(2^) 
,r(i + 2/i) 



[2isin7r(/i - u>) + e^"^ 1 "^*] 



T(1 + 2 M ) 



[ e -«(M+<s)2isin7r2/i + e ™(-^+a)£*] 



r(l + 2M) p 7ri(-^+gi) £*\ 



Miv = 



r(2m) 



r(2 Wl ) 



'(1 + e ~™ 2ui [2^sin7r(/z + cj) + e -"-*(/*+u0^*] g-TriC/H-w)^ _ £*) 



3 7r!(2/j-2ui) 



2?ri 



,r(i + 2/i)r(i - 2wi) 



2ni 

r(-2 / /)r(2w 1 ) 

e «(2^-2wi) 2costt(2/x + 2wi) - e^ 2 ^ 2 " 1 ) 



> 

> 

a 

> 

a 
3 



H 

a 
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Proof. The set of Pyi parameters p.4j) in our application imply that the limiting 
Py parameters are 9q = — 2u>i, 6*o v = H + &, #iv = —^ — w, #oov = 2/i — 2wi, which 
are consistent with those of (|5.8[) . Using the structure of the Pyi monodromy 
matrices (J37TJ) and (|3.8p wc compute that ex — — ^oov ^6 and /? — — $6, which fixes 
the Stokes multipliers as given by (I6.19p . We solve the system (|6.17|) which yields 
the solution 



for an arbitrary non-zero K\ \. Finally we employ these findings in (|6.16p and 
after lengthy calculations arrive at explicit forms for the monodromy matrices 




(6.23) 



A" 



2ig7ri(0ooV+0 6 ) 

Si 



(lG.20lG.21l6.22h . 



□ 



However the explicit parameterisation of the alternative Pv monodromy matrices (|4.14I4.15|) allows us to compute these directly 
from the data deduced from the expansion of the T-function about x = 0, as given by Theorems 15.11 and 15.21 

1 [5] In terms of the other set of monodromy data {Mov, Miv, M^y} we have 

e 7ri(-3M+o>)(i_£*) _ £i^+|/f) e ^(-^) [2isin7r2^ + e-" 2 ^C] 

[2i sin7r(/i - u>) + e ^(-M+ffl)£*] 



Corollary 6.2. 



(6.24) 



M, 



ov 



r(2^x) 



and 
(6.25) 

and 
(6.26) 



\r(i + 2/i) e 
/ 



2cos7r(^ + w) - e "(- 3 ^)(l - £*) 



M 



IV 



r(i + i( _ u+uj)£:t 



M 



ooV 



'2costt(2/i + 2wi) _ e -^(2M-2wi) 



2ttz 

r(-2/i)r(2wi) 

27TZ _-?ri(2/i— 2wi) -7ri(2^-2wi) 

r(i + 2/i)r(i - 2wi) 

Proof. The first two matrices are computed from (|4.14p and (|4.15p after noting that a comparison of the Stokes multipliers (14. 1 7[) 
and (|6.19p obliges us to set r = — 2/x. The third matrix is computed using (|4.7p . As a check one can verify directly that both sets 
of monodromy matrices (|6. 2016. 2116. 22|) and (|6. 2416. 2516. 26[) are related by the transformations of (|4.9p . □ 
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In the discussion of the asymptotics of the Py cr-function as t — > ±00 given in 
[2] the following parameter 

(6.27) 0o = ^ln{f[l - e* 1 (- a H-a-i)(i _£*)]) 

was found to be crucial. We note that when £* = 1 both /3q and the upper left 
element of Mq both vanish. A comparison of our monodromy matrices, evaluated 
at £* = 1, with those discussed in p] where a solution to one of the connection prob- 
lems was reported reveals that our case is precisely the case of the lower truncated 
solution. In terms of the parameters used in that work we have 

2(2^i-2wi) -7ri(^-i>) ^(l + 2/x) 



(6.28) iu = 2 2 W- 2 " 1 >e 



r(2wi) ' 



(6.29) iv = J~e^ +Ul) cos7r(/i-u>). 

This means that the cr-function has the following asymptotic expansion 

(6.30) C(s) ~ ( (s)~ie m ^ + ^cosn(ti-Lu)\M-e- ls / 2 

s — >— ioo V 47T 

in the sector — 7r < arg(s) < 0. Here (o(s) is the formal, algebraic expansion 

s 2 1 

(6.31) Co(s)~ tt + (m-o*w 2 )s 

s|^oo 16 2 

+ 4fi 2 - 2^iw 2 + + lu' 2 - - - 2iw 2 (4^ 2 - 4o> 2 )s~ 1 
+ [l6/z 4 - 8(4(w 2 + wfj + IV 2 - 16w 2 w| + (4cj 2 - l)(4w 2 - 4a; 2 . - 1)] S - 2 +0(s" 3 ). 
However this is not the physically interesting asymptotic expansion of s — > +ioo. 

7. Conclusions 

An especially interesting problem that remains outstanding is the connection 
formula relating the x — > behaviour of the Py r-function arising in the bulk scaling 
of the spectrum singularity ensemble A(x), as given in (|5.7jl . to the x — > +ioo 
behaviour. The bulk scaling limit of the partition function for the Dyson CUE 
(|3.ip is the generating function for the gap probability, which has the Fredholm 
determinant formula 

(7.1) E((-t, t);0 = det(I - £K\ LH _ t>t) ), 

where the integral operator IK has a kernel with the sine kernel form 

sin7r(x — y) 

7.2) K(x,y) = —r± 

n{x - y) 

This is a special case of our spectral average A(x) and the precise relation is 
(7.3) E{(-t, t);0 = A^t; M = 0, wi = 0, u 2 = 0; £). 
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The large t — > oo asymptotic expansion of the gap probability E((— t, t);£ = 1) is 
known to be 

(7.4) E((-t,t);£ = 1) ~ e 3C'(-i)+Aiog2 t -i/4 e -i t 2 +o(i) 5 
or alternatively as 

(7.5) tj t logE((-t,t);Z = l)~-t 2 
For < £ < 1 this result becomes 

(7.6) i^log£?((-t,t);0~-log(l 

at 7r 

This asymptotic problem has been the subject of intense and continuing study 
PHni|31IH51IHl[niinilZl- The question we pose then is - what is the generalisation 
of this asymptotic result valid for generic Wxj^jA* an d that is uniformly valid for 
|1 — £| < 6 with some finite 5 > 0? 

Our results do not strictly apply for the situation of moments or singularities 
that are negative 5R(^i) < —1/2 and \t\ — > 1, however the present results 

may carry over to this case. This case is quite relevant in the context of studies 
concerning the averages of the characteristic polynomial in the CUE with negative 
integer powers [12] where such averages are employed to make conjectures concerning 
the averages of powers of the zeta or L- functions on their critical lines. 
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